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ABSTRACT 

We investigate the magnetic geometry of the active G8 dwarf £ Bootis A, from spec- 
tropolarimetric observations obtained in 2003 with the MuSiCoS echelle spectropo- 
larimeter at the Telescope Bernard Lyot (Observatoire du Pic du Midi, France). We 
repeatedly detect a photospheric magnetic field, with periodic variations consistent 
with rotational modulation. Circularly polarized (Stokes V) line profiles present a sys- 
tematic asymmetry, showing up as an excess in amplitude and area of the blue lobe of 
the profiles. A direct modeling of Stokes V profiles suggests that the global magnetic 
field is composed of two main components, with an inclined dipole and a large-scale 
toroidal field. We derive a dipole intensity of about 40 G, with an inclination of 35° of 
the dipole with respect to the rotation axis. The toroidal field strength is of order of 
f 20 G. A noticeable evolution of the field geometry is observed over the 40 nights of 
our observing window and results in an increase of the field strength and of the dipole 
inclination. 

This study is the first step of a long-term monitoring of £ Bootis A and other active 
solar-type stars, with the aim to investigate secular fluctuations of stellar magnetic 
geometries induced by activity cycles. 

Key words: Line : polarization - Stars : rotation - activity - magnetic fields - Star : 
individual : HD 131156 



1 INTRODUCTION 

Fluctuations in the magnetic activity of the Sun are domi- 
nated by the well-known solar cycle, historically discovered 
as a pall yr period of the sunspot number (Schwabe 1843), 
which itself was later linked to a polarity reversal in bipo- 
lar magnetic regions between both hemispheres of the Sun. 
The global activity cycle of the Sun, involving two consec- 
utive polarity switches, is therefore of period «22 yr (Hale 
& Nicholson 1925) . Efforts to develop a theoretical explana- 
tion of such a cyclical activity have given rise to the so-called 
solar dynamo theories (see e.g. Ossendrijver 2003 for a re- 
view) . 

In solar-type stars, long-term monitoring of chromo- 
spheric emission suggests that the existence of activity cy- 
cles depends a lot on the evolutionary stage and rotation 



rate (Baliunas et al. 1995). It is suggested in particular that 
a cyclical activity is common among slowly-rotating, low ac- 
tivity main sequence stars similar to the Sun. By contrast, 
younger stars (rotating faster and possessing a much higher 
level of activity) generally display more erratic fluctuations 
of activity level. If one is to understand the dynamo pro- 
cesses at work in solar-type stars, and how various stellar 
parameters may affect the onset of different types of dy- 
namos, it is necessary to obtain direct information about 
the geometry of stellar magnetic fields and about secular 
fluctuations of this geometry. 

For rapidly rotating stars, the surface distribution of the 
magnetic field vector is now routinely reconstructed with the 
technique of Zeeman-Doppler Imaging (ZDI hereafter, Do- 
nati & Brown 1997, Donati et al. 2003) . Mapping techniques 
are however much more limited for slow rotators, because in 
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their case all magnetic regions of the visible hemisphere of 
the star produce Zeeman signatures of similar radial veloc- 
ities. Because of the usually complex field topology of cool 
active stars, spectroscopically superimposed Zeeman signa- 
tures of mixed polarities are observed at any time, resulting 
in an important cancellation of the polarized signal. This 
situation is occurring for slow rotators of course, but also 
for low-mass stars with small radii and for stars with low 
inclination angles. For such objects, most individual active 
regions remain unresolved. The only measurable quantity 
is the global field of the star, i.e. the magnetic field that 
remains after integration over the visible stellar hemisphere. 

In the case of the Sun, the large-scale structure of the 
photospheric field can be determined using magnetograph 
data (e.g. Stenflo 1991). Outside the belts of active regions, 
and in particular around the poles, a large-scale component 
of the surface field is observed (sometimes called the back- 
ground field of the Sun). This component was first investi- 
gated by Babcock & Babcock (1955), who report that its 
geometry is mostly dipolar and close to axi-symmetry. The 
solar background field is modulated by the Hale cycle, in 
anti-phase to the sun-spot cycle (its polarity reversal occur- 
ring at solar maximum). The strength of the field at the 
pole during solar minimum is of the order of 5 G (Smith & 
Balogh 1995). Since its cyclical evolution is controlled by the 
large-scale dynamo of the Sun, investigating the geometry 
and secular evolution of global fields on a sample of slowly- 
rotating late-type stars may also bring precious information 
about the dynamo processes at the root of their magnetic 
activity. 

As a first attempt in this direction, we concentrate on 
the solar-type star £ Bootis A. From photospheric and chro- 
mospheric activity tracers, this object can be considered 
as a very active star, with irregular fluctuations of activity 
(Toner & Gray 1988, Baliunas et al. 1995, Gray et al. 1996). 
Its surface magnetic field was first detected by Robinson 
et al. (1980). The rotational modulation of the large-scale 
field was then investigated by Plachinda & Tarasova (2000) 
from a time series of longitudinal field measurements col- 
lected over some two decades. We propose here a different 
and potentially more powerful approach, making use of a 
time-series of high signal-to-noise spectropolarimetric data, 
and based on a direct modeling of the polarized line profiles, 
from which we derive tight constraints on the field geometry. 

The article is divided as follows. We first discuss the 
fundamental parameters of £ Bootis A. We then present our 
set of spectropolarimetric data and investigate the large- 
scale geometry of the photospheric field. We discuss our re- 
sults in relation with previous works on magnetically active 
stars and propose some directions in which this study may 
be extended in a near future. 



2 FUNDAMENTAL PARAMETERS 

£ Bootis A is a bright member (m v = 4.7) of one of the 
nearest visual binary systems (parallax it — 0.1491 trans- 
lating into a distance of 7 pc, Gliese & Jahreiss 1991). The 
system has been regularly observed for decades and is now 
known to follow a 151 yr orbital period. The masses of 
£ Boo A and B are estimated to be 0.85 and 0.72 Mq re- 
spectively and the inclination angle of the orbital plane is 



140° (Wielen 1962). The respective spectral types of both 
components is G8V and K5V (Abt 1981). The mean effec- 
tive temperature of the primary component (T e ff ~ 5550 K) 
was determined by Gray (1994), using line-depth ratios cor- 
rected from metalicity dependence. Noticeable rotational 
and secular fluctuations of T c s, presumably produced by 
magnetic activity, are also reported by Gray et al. (1996) 
and reach a level of 12 K. The metalicity of the primary 
dwarf is [Fe/H] = -0.20 ± 0.08 (Cayrel de Strobel et al. 
1992) and its luminosity is Log(L/L Q ) = -0.26 ±0.03 (Fer- 
nandes et al. 1998). Using stellar evolutionary models for 
both stars of the system, Fernandes et al. (1998) derive dif- 
ferent additional stellar parameters. They propose an helium 
abundance Y = 0.26 ± 0.02 and a mixing length parameter 
Q.MLT = 1.7±0.03. Their models are consistent with a young 
age (2 ± 2 Gyr), which is also suggested by the high chro- 
mospheric emission (Baliunas et al. 1995) and high lithium 
abundance (Herbig 1965). 

The rotation period has been estimated in several stud- 
ies, making use of different rotationally modulated tracers 
(Ca II emission by Noyes et al. 1984, line asymmetry and 
line ratios by Toner & Gray 1988, longitudinal magnetic field 
by Plachinda & Tarasova 2000). All periods are in the 6.1- 
6.5 d range, but the tight error bars on each individual es- 
timate reveal an apparent discrepancy. The shortest period 
is derived from the integrated longitudinal magnetic field 
(P~6.1455 d). The longest one is obtained from variations 
in line bisectors (P«6.43 d). We use this last period, care- 
fully derived by Toner & Gray (1988) from a high-quality 
data set, as our reference to determine the rotational phases 
of our observations. We will show in Sect. 14. ll that this esti- 
mate is also consistent with the rotation period derived from 
our own observations. We take Julian date 2,452,817.41 (the 
date of our first observation) as phase zero. 

From the luminosity estimated by Fernandes et al. 
(1998), we deduce a stellar radius R = 0.8 ± 0.03 R Q . 
Combining this value with the projected rotational velocity 
v sin l= 3±0.4kms _1 (Gray 1984), we obtain an inclination 
angle i = 28 ± 5°, which gives also the equatorial velocity 
w e = 6.4 ± 0.3 kms^ 1 . In previous studies (Toner & Gray 
1988 and after them Plachinda & Tarasova 2000), the ro- 
tation axis was assumed to be perpendicular to the orbital 
plane of the binary system, i.e. i was chosen close to 40°. 
However, considering the very long orbital period of the sys- 
tem, the hypothesis of a tidally-induced perpendicularity of 
the rotation axis with respect to the orbital plane is not ex- 
pected (Zahn 1977). We therefore prefer our own estimate 
of i, which is adopted in the rest of the paper. 



3 INSTRUMENTAL SETUP AND 
OBSERVATIONS 

Observations were carried out at the Telescope Bernard 
Lyot, using the MuSiCoS echelle spectrograph (Baudrand 
& Bohm 1992) associated to its Cassegrain-mounted polari- 
metric module (Donati et al. 1999). 

In polarimetric mode, MuSiCoS spectra cover the 450- 
660 nm range, with a resolving power close to 30,000. 
The analysis of circular polarization is performed with a 
quarter-wave plate (transforming circular polarization into 
linear polarization), followed by a Savart plate which splits 
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Table 1. Journal of observations. From left to right, we list the date of observations, the Julian date, the UT time, the S/N of raw 
Stokes V spectra and of LSD mean profiles, the rotation cycle (according to a 6.43 d period and taking our first observation as phase 
zero), the longitudinal field B[ and the uncertainty as on the longitudinal field measurement. 
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the incident light into two beams, respectively containing 
light linearly polarized perpendicular/parallel to the axis of 
the beam-splitter. Each circularly polarized spectrum (also 
called Stokes V spectrum hereafter) is obtained from a se- 
ries of four sub-exposures taken with the quarter-wave plate 
oriented at azimuth ±45° with respect to the optical axis of 
the Savart plate. Following this procedure (detailed by Do- 
nati et al. 1999), both beams exchange their optical paths 
and positions on the detector, ensuring a removal of most 
spurious signatures to less than 5 x 10 _5 / c rms, where 7 C 
denotes the continuum intensity. The intensity (Stokes I) 
spectrum is simultaneously produced by a simple addition 
of all four frames. Using a third combination of the series of 
sub-exposures, a "null spectrum" is obtained (iV hereafter, 
Donati et al. 1997), which does not contain any polarized 
signal. 

The data set was collected during the summer of 2003, 
with 28 Stokes V spectra gathered from June 26 to August 
04 (Table 0. The data reduction was performed using the 
ESpRIT package (Donati et al. 1997). In order to improve 
the wavelength calibration usually obtained for MuSiCoS 
spectra with reference to a Thorium- Argon spectrum, an 
additional calibration was performed, making use of atmo- 
spheric lines over-imposed on the stellar spectra. Using a 
line mask containing some 470 telluric lines, we thus reach a 
radial velocity accuracy of order of 300 m.s -1 (Petit et al. 
2004a). The exposure time of each individual sub-exposure 
is 500 s, leading to a signal-to- noise ratio (S/N) of polarized 



spectra of the order of 550 in optimal weather conditions. In 
order to decrease further the noise level, the Least-Squares 
Deconvolution technique (LSD, Donati et al. 1997) was used 
to generate high S/N line profiles by simultaneously extract- 
ing the information from all photospheric lines available in 
the echellograms. Using a line mask calculated for a stellar 
atmosphere of spectral type G7, about 2,700 lines are avail- 
able in the reduced spectra, allowing us to produce mean 
line profiles with S/N close to 20,000 (i.e. a typical S/N 
multiplex gain of 35 with respect to the peak S/N of the 
raw spectra). 

The series of Stokes V profiles is plotted in Fig.Q in the 
rest-frame of the star, i.e. after correction of the heliocentric 
velocity (varying around an average value of -24 km s~ ) and 
of the mean stellar velocity (equal to 1.9 kms -1 at the date 
of observations) . Zeeman signatures are repeatedly detected 
inside line profiles, always with a classical double-peaked 
shape appearing in the ±10 kms~ velocity range. The am- 
plitude of signatures varies between 0.02% and 0.07% of the 
(non polarized) continuum level. The highest signal mea- 
sured (on Jul. 18) has a significance of about 10<r, where a 
represents the uncertainty on the measurement. No spurious 
signature is detected in the N profiles (Fig. Q right panel), 
confirming the robustness of the field measurements. We 
note also that spectra of magnetic and non-magnetic stan- 
dard stars were obtained during the same observing run, 
verifying the nominal functioning of MuSiCoS. The global 
shape of Zeeman signatures is mostly stable with time (no 
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Figure 1. Rotational modulation of line profiles. Left : Stokes V profiles calculated with a G7 line-mask. Right : "Null" profiles. 
White/black saturation levels correspond to ±3 X 10~ 4 I C . 



Mean Stokes I profile Mean Stokes V profile 
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Figure 2. Mean line profiles calculated from an average over all available observations. Left : Stokes I profile. Right : Stokes V profile 
(full line) and N profile (dashes). The dotted line indicates the zero level. 



polarity switch observed during the rotation, for instance), 
with a trend for their amplitudes to be larger around phase 
0.5. 

A systematic asymmetry is observed between both lobes 
of the Stokes V profiles. If we call respectively a b and a r the 
amplitude of the blue and red lobe of the averaged profile 
(plotted in Fig. |2j, we derive a relative asymmetry 8a = 
(at — a r )/(a,b+a r ) = 0.19±0.05. A difference is also visible in 
the areas of the lobes, with a relative asymmetry S A = (At — 
A r )/(A b + A r ) = 0.22 ±0.04 (where A b and A r are the areas 
of the blue and red lobe, respectively). No counterpart is 



seen in the mean N profile, suggesting that such asymmetry 
is not of instrumental origin. We note that no obvious phase 
dependence of 5a and 5 A is observed, owing to the limited 
S/N of the data set. 

The averaged Stokes I LSD profile is plotted in Fig. [5] 
The line asymmetry investigated by Toner & Gray (1988) 
and Gray et al. (1996) using line bisectors is also observed 
in our data set, despite much lower spectral resolution. Bi- 
sectors are plotted in Fig. El f° r LSD profiles averaged over 
different bins in rotational phase. Each bin covers 20% of 
the stellar rotation. We note that bisectors display a sys- 
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Figure 3. Bisectors of Stokes I LSD profiles averaged over differ- 
ent bins in rotational phase. Each bin covers 20% of the rotation 
cycle, from phases 0.0-0.2 to phases 0.8-1.0. 



tematic shift toward the red near the continuum, which is 
consistent with previous findings of Toner & Gray (1988). 
If we define the velocity span of a bisector as the difference 
in radial velocity between a point near the top of the profile 
(I/I c = 0.95) and a point close to the bottom (I/I c = 0.8), 
we find that the velocity span is marginally higher around 
phase 0.6 (with a span of about 170 ms -1 , vs. 110 ms -1 
around phase 0.9). All other types of variability observed 
in Stokes I profiles are close in amplitude to the expected 
instrumental errors. 



4 GEOMETRY OF THE LARGE-SCALE FIELD 

4.1 Longitudinal field measurements 

For every Stokes V and I profiles of the data set, we calculate 
the corresponding longitudinal magnetic field Bi using the 
following expression (Donati et al. 1997, Wade et al. 2000) : 



Bi = -2.14 x 10 1 



J vV(v)dv 



Xozc f (I c — I(v))dv 



(1) 



where v is the radial velocity, Ao the mean wavelength of the 
line mask used to compute the LSD profiles (521 nm in our 
case) and z the mean Lande factor of the line-list, equal to 
1.21. The integration window covers a ±22 kms" 1 velocity 
range around the line centroid. The derived values of Bi 
are listed in Tabled The longitudinal field varies between 
-2 G and 18 G, with uncertainties ranging from 3 to 10 G. 
Significant fluctuations show up around the mean value of 
8 G. 

The Bi time-series was used for a period search. The 
rotation period we obtain is 6.36 ± 0.1 d, in agreement with 
the period derived by Toner & Gray (1988). The fact that no 
more than 6 consecutive periods are covered by our data-set 
explains the relatively large error bar we get here. 

We plot in Fig. [I] the rotational modulation of Bi, ac- 
cording to the 6.43 d rotation period adopted in the paper. 
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Figure 4. Longitudinal field as a function of the rotational phase 
for data collected until July 17 (filled squares) and after this date 
(open squares). The dashed curve represents a least-square fit of 
Preston's equation for the whole data set. The dash-dotted line 
illustrates a null magnetic field. 



In the case of a tilted dipole, the rotational modulation of 
Bi is described by the equation of Preston (1967) : 



B, =B 



15 + u 



p 20(3 



(cos j3 cosi + sin j3 sini cos[2-7r(0 — 0o)])(2) 



where B p is the field strength at the magnetic pole, (f>o the 
rotational phase at which the observer faces the positive 
magnetic pole and j3 the inclination angle between the rota- 
tion and magnetic axis, it is the limb-darkening coefficient, 
that we take equal to 0.66 according to Wade & Rucinski 
(1985). The inclination i is set to 28°. Making an iterative 
adjustment of Preston's relation, we obtain B p = 41 ± 5 G, 
(3 = 33 ± 6°and (p = 0.59 ± 0.04, with a reduced X 2 (X? 
from now on) equal to 1.2 (vs. x?=4.5 when assuming the 
field to be constant and equal to its mean value 8 G during 
the whole rotation). 



4.2 Direct fitting of Stokes V profiles 

Rather than restricting the modeling of the field geometry 
to the first moment of Stokes V profiles, we conduct in this 
section a direct modeling of the profiles. 

With this aim, we first check that the adopted val- 
ues of stellar parameters enable a correct description of the 
Stokes I LSD profiles. The line model we adopt is similar 
to that described by Donati & Brown (1997), except that 
the smearing induced by macro-turbulence is not assumed 
to be Gaussian in our case. A good fitting of stellar line pro- 
files is obtained with vs'mi—3 krns" 1 and macro-turbulence 
Cut = 5 kms -1 , taken from Gray (1984). 

As a second step, we generate a synthetic time-series 
of Stokes V LSD profiles from an artificial star possessing 
a centered dipolar field, with rotational phases of the syn- 
thetic profiles taken equal to that of our observations. The 
model used to calculate polarized profiles is again that of 
Donati & Brown (1997) where, for each pixel of the stellar 
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surface, the local Stokes V profile is assumed to be propor- 
tional to the derivative of the local Stokes I profile. We then 
tune the values of B p , (f>o and /3. Using a total of about 1,700 
different values of the three field parameters, we generate as 
many synthetic data sets that we compare to the observa- 
tions. A x 2 is calculated for each model, so that we end up 
with a three-dimensional \ 2 landscape. The location of the 
minimum x 2 gives the most likely value of the parameter 
triplet. 

Around the \ 2 minimum, error bars on the parameters 
(considered separately) correspond to an increase of x 2 of 
unity (Press et al. 1992). Each error bar is multiplied by 
(Xmin) 1 ^ 2 ! m order to take partially into account potential 
systematic errors of the method (e.g. departure of the ob- 
served photospheric field from an inclined dipole). 



4.3 Global dipole 

Following this procedure, we derive B p = 35 ± 2 G, (j>o = 
0.58±0.04 and f3 = 22±5° (Tab.gJ. The minimum x? is not 
better than 2.3, suggesting that our model provides nothing 
more than a rough approximation of the profiles. In Figure 
|S] we represent 2D-maps that we extract from the 3D-x?- 
landscape by fixing one of the three free parameters to its 
best value. By doing so, we obtain a convenient visual con- 
trol of the fitting procedure, showing that the x? minimum 
is unique, with well-defined ellipses of iso-x? surrounding the 
minimum. The synthetic profiles (averaged over the stellar 
rotation) are compared to the observations in Fig. |H] 

As a second test of the robustness of our method, we run 
again the x? minimization, but this time using the N profiles 
(following Wade et al. 2005). As expected, the derived field 
strength B p is now compatible with zero (Tab. |2J and all 
values of (3 and tpo give similar x? f° r a fixed value of B p . 
Another major difference with the fitting of Stokes V profiles 
is that the minimum x? is now close to unity. The fact it is 
actually less than one (0.92 at minimum) suggests that error 
bars calculated for MuSiCoS spectra are over-estimated by 
about 4% (as already pointed out by Wade et al. 2000). 
It tells also that the relatively high Xr obtained with the 
Stokes V data set is probably due to intrinsic limitations in 
our magnetic model. 

At this stage, we note that the data fitting is not im- 
proved by going further in the multi-pole development of 
the field. If we add a quadrupolar field component (as- 
sumed to be aligned on the original dipole), the best fit is 
consistent with a null quadrupole (with a derived strength 
B q = — 5 ± 5 G). The x? is not significantly decreased by 
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the addition of the quadrupole. 



Figure 6. Synthetic Stokes V profiles (thick lines) compared to 
observed profiles (thin lines), in the case of a simple dipolar fit. 
Rotation cycles are indicated right to each profile. Successive pro- 
files are shifted vertically, for display purpose. The dotted lines 
illustrate the zero level of each profile. 



close to 2.3, with variations not exceeding 2% for all three 
values of i. 

A x^Tandscape is also created for a rotation period of 
6.14 d. Again, the Xr is varying by less than about 1% with 
respect to the original period of 6.43 d. New values of B p 
and P keep within the error bars of the original values. Be- 
cause the rotational phases are calculated with the new value 
of the period, <f>o changes significantly and is now equal to 
0.74 ±0.05. 

From both tests, it is clear that the limited quality in 
data fitting we obtain with a purely dipolar field is unlikely 
to come from a wrong estimate of stellar parameters. 



4.4 Impact of stellar parameters 

In this section, we investigate whether a possible error in 
the determination of stellar parameters may contribute to 
the relatively high level of the x? we obtain with the fit of 
a dipolar field. 

We first vary the inclination i. The main outcome of 
this test is to show that the magnetic model does not de- 
pend critically on this parameter. The field geometries we 
calculate for i = 20° and i = 40° are in good agreement 
with that we get for i = 28° (Tab.|2J. The best x? is still 



5 LARGE-SCALE TOROIDAL FIELD 

Since a potential error in fundamental stellar parameters 
does not seem to be at the origin of the poor data adjust- 
ment, we propose in this section and in Sect. |H|a series of 
tests, in order to understand better how various effects can 
lead to the observed limitations of the dipole model. Differ- 
ent options are therefore explored to improve the modeling 
of the magnetic geometry, and for each test we discuss the 
impact observed on the values of the field parameters, as 
well as on the value of the x?- 
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Figure 5. maps obtained from slices in the three-dimensional field parameter space. B p is the magnetic field strength at the magnetic 
pole. (3 is the inclination of the dipole with respect to the rotation axis. </>o is the rotational phase of maximum field strength. 



Table 2. Estimated magnetic parameters for the various models described in the paper. Values without error bars are taken constant 
during the x 2 minimization process. The best-fitting models are represented as bold rows. 



Bp <j> B t Xr Notcs 

(G) (°) (G) 



35 ± 2 


0.58 ±0.04 


22 ± 5 





2.35 


dipole only 


0±2 









0.92 


dipole, N profiles 


33 ± 2 


0.59 ± 0.05 


22 ± 7 





2.37 


dipole, inclination i = 20° 


39 ± 2 


0.59 ± 0.05 


18 ±7 





2.30 


dipole, inclination i = 40° 


34 ±2 


0.74 ±0.05 


21 ±7 





2.32 


dipole, rotation period of 6.14 d 


33 ± 2 


0.48 ± 0.07 


23 ± 7 





2.26 


dipole, odd-numbered profiles 


31 ± 1 


0.45 ± 0.05 


22 ± 6 





2.29 


dipole, even-numbered profiles 


27 ±2 


0.64 ±0.13 


9 ± 7 





1.92 


dipole, first half of the data set 


40 ±2 


0.55 ± 0.03 


28 ±4 





2.53 


dipole, second half of the data set 


39 ±2 


0.58 


35 ±3 


-148 ±9 


1.30 


dipole + toroidal field 


0±2 


0.58 




0±8 


0.92 


dipole + toroidal field, N profiles 


30 ±2 


0.58 


24 ±4 


-124 ±8 


1.22 


dipole + toroidal field, first half 


43 ± 2 


0.58 


35 ± 3 


-160 ± 8 


1.33 


dipole + toroidal field, second half 


35 ± 2 


0.58 


31 ±4 


-137 ±8 


1.36 


dipole + toroidal field, odd-numbered prof. 


37 ±2 


0.58 


35 ±4 


-149 ±8 


1.32 


dipole + toroidal field, even-numbered prof. 


37 ±2 


0.58 ± 0.04 


23 ± 5 





2.31 


dipole, 2-component atmosphere 


37 ±2 


0.57 ±0.05 


22 ± 6 





2.60 


dipole, rising flow only 


401 ± 26 


0.63 ± 0.03 


36 ±4 





2.31 


dipole, sinking flow only 


38 ± 2 


0.58 ± 0.03 


23 ±4 





1.87 


dipole, asymmetric synthetic profiles 


40 ±2 


0.58 


33 ±3 


-112± 8 


1.23 


dipole + toroidal field, asym. profiles 



5.1 Large-scale toroidal field 

Since a simple dipolar geometry of the surface field is not 
able to produce a satisfactory data fitting, we propose to 
modify our magnetic model by adding a new component to 
the large-scale photospheric field. We choose to impose the 
presence of an azimuthal field component (in addition to 
the inclined dipole). This addition is inspired by a result of 
ZDI, revealing the existence of large-scale azimuthal fields 
on most extremely active stars observed to date (e.g. Donati 
et al. 2003, Petit et al. 2004b). 

As a very rough assumption, we take the strength of 
this component to be constant over a ring encircling the 
pole between latitudes 30° and 60° , which is somewhat rem- 
iniscent of the large-scale structure of this component on 



several rapidly-rotating stars observed so far. Because this 
ring possesses an axi-symmetric geometry, we assume that 
the determination of (j>o will not be affected by its pres- 
ence. We therefore keep 0o equal to 0.58 and adjust again 
three parameters simultaneously : B p , [3 and the toroidal 
field strength B t . 

The best fit is obtained for B t = -148 ± 9 G (where 
B t < corresponds to a clockwise field). The value of B p is 
slightly affected by the presence of the azimuthal field (39 ± 
2 G). P is also slightly increased to 35 ±3°. Slices taken from 
the 3D-x 2 -landscape are represented in Fig. |7| We obtain 
here a very significant improvement in the quality of the fit, 
with Xr now equal to 1.3. As previously done with the dipole, 
we check that the same \ 2 minimization, applied to the N 
profiles, results in a toroidal component consistent with zero 
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(Tab.© and a Xr close to unity. The value of the azimuthal 
field depends of course on the latitudinal extent of the ring 
on the synthetic star (mainly as a result of the different 
filling factor). Assuming for instance that the ring fills the 
whole hemisphere of the star reduces the field strength to 
B t — — 81 ± 5 G, but no significant impact is observed on 
the other field parameters, as well as on the minimum x?- 



5.2 Short-term evolution of the field geometry 

Despite very significant improvement with respect to a 
purely dipolar field model, a x? of unity is still not reached 
with the additional assumption of a toroidal field compo- 
nent. In this section, we investigate to what extent an intrin- 
sic evolution of the field during data collection may produce 
such limitation. 

In order to investigate potential short-term changes in 
the field geometry, we generate a couple of sub-sets from our 
original data set, the first one containing observations of the 
first half of the run (until July 17), while the second one 
contains observations gathered after this date. Both sub- 
sets present a correct phase sampling, as can be checked 
from Tab. [U 

Assuming first a purely dipolar geometry, we determine 
the field parameters for each new data set. From the first 
one, we derive a field strength B p = 27 ± 2 G and an incli- 
nation of the dipole axis (3 — 9 ± 7° . The parameters derived 
from the second one are B p = 40 ± 2 G and (3 = 28 ± 4°. 
Both sub-sets give consistent values of (f>o- The minimum x? 
is better for the first set (1.9 versus 2.5). We suggest that 
this apparent evolution is real and results from an intrinsic 
evolution of the stellar magnetic field over the 40 consecu- 
tive days of our observing window. A direct look at the Bi 
measurements confirms this evolution. As an example, Bi is 
equal to Bi = 18 ±4 G and 15 ±5 G respectively at rotation 
cycles 4.811 and 5.747 (July 27 and Aug. 02), while the lon- 
gitudinal field does not exceed — 2±5 G at phase 1.707 (July 
07). We note that changes in the field geometry during data 
collection are not responsible for the high Xi 01 the dipole 



model, since the Xi level is not significantly decreased by 
the splitting of the original data set. 

As for the global dipole, we check if an evolution is 
also observed for the toroidal component. Using again the 
same couple of subsets, we calculate two Xr-landscapes and 
optimize B p , /3 and -Bt. We report an increase of Bt, from 
-124 ± 8 G to -160 ± 8 G. The increase of B p and (3 is 
also observed here. Again, the splitting of the data set is 
not sufficient to provide a Xi of unity (x? = l-22 and 1.33 for 
the first and second sub-sets, respectively). 

In order to double-check the reality of the field evolu- 
tion, we now split the original data set into a new pair of 
sub-sets, the first one containing odd numbered profiles and 
the second one built up with even-numbered profiles. The 
dipole and toroidal field parameters obtained in both mod- 
els are in close agreement with the values derived from the 
full data set (Tab.©. 



6 EFFECT OF SURFACE VELOCITY FIELDS 

In Fig. [HI we plot the derivative of the averaged Stokes I line 
profile, together with the mean Stokes V profile. The Zee- 
man signature presents a clear velocity shift with respect to 
the line centroid, showing up as a 2.5 kms -1 red-shift of 
the velocity at which the profile changes sign ( called zero- 
crossing velocity from now on). This shift is well above the 
300 ms _1 radial velocity accuracy of the spectrograph. It 
reaches a maximum of about 5 kms -1 at rotation cycles 
1.239, 1.707 and 6.056, but no evidence of a smooth phase 
dependence is observed. In this context, a simple magnetic 
dipole (producing Zeeman signatures roughly proportional 
to dl /d\) cannot reproduce the observed shape of Stokes V 
profiles. A large-scale azimuthal field component mostly af- 
fects circularly polarized profiles when seen at intermediate 
limb angle. With an axi-symmetric ring of azimuthal field, 
two main regions of different radial velocities (on opposite 
sides of the meridian facing the observer) contribute to form 
two Zeeman signatures of opposite signs. The whole ring 
therefore produces a Zeeman signature resembling d 2 I/d\ 2 
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Figure 8. Derivative of the averaged Stokes I profile (full line), 
compared to the averaged Stokes V profile (dashes). 

and is able to compensate for the apparent velocity shift of 
Stokes V profiles. 

In this section, we investigate to what extent a system- 
atic velocity of magnetic regions (presumably related to sur- 
face convective flows) may also contribute to the observed 
red-shift of circularly polarized profiles. 

6.1 Two-component model of convective motions 

We adopt a modeling of stellar granulation inspired from 
the work of Gray & Toner (1985). The surface convective 
flows are described in a two-component model, with hot 
(rising) and cool (sinking) material. Each component is as- 
sumed to produce a Gaussian line profile of dispersion £, cen- 
tered around radial velocity « C ooi and iw for the down-flow 
and up-flow respectively. The line-depth ratio between both 
Gaussian components is quoted R c /^ thereafter. The granu- 
lation parameters we adopt are those Toner & Gray (1988) 
derived specifically for £ Bootis A, with £ = 5 kms" 1 , 
"cool — "hot = 5.4 kms -1 and R c /h = 0.1. The adopted 
value Vhot = —0.4 kms -1 is the only parameter we modify 
from the value proposed in the model of Toner & Gray (in 
which "hot — —1.4 kms" 1 ). A slight red-shift of 1 kms" 1 
with respect to the original model is imposed, in order to 
take into account that the predominant hot component pos- 
sesses a radial velocity close to zero in our data set (as a 
result of the correction of the mean radial velocity of the 
star). The same red-shift is imposed to " CO oi, in order to keep 
"cool — "hot equal to the value of the original model (bearing 
in mind that this parameter, along with R c /h, controls the 
final shape of the line profile). 

The main effect of this new line model is to introduce 
an asymmetry in the shape of Stokes I profiles. The level of 
asymmetry is however small and does not produce a signif- 
icant impact in Stokes I fitting. For Stokes V, if we assume 
at first that the magnetic dipole is the same in granules and 
in inter-granular lanes (i.e. identical Zeeman signatures are 
formed in both components of the atmosphere), the data ad- 
justment leads to a Xr very similar to that obtained in Sect. 
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Figure 9. Averaged Stokes V profiles calculated with blue and 
red line-masks (full line and dashes, respectively). 

14.31 (see Tab.|2J. The best-fitting dipole parameters (Bp, f3, 
4>o) remain also very close to the original values of Sect. 14.31 

As a second step, we run again the \ 2 minimization 
procedure, but this time we keep only one of the two com- 
ponents of the convective atmosphere. By doing so, we as- 
sume that magnetic regions are concentrated either in the 
rising or sinking material alone. Considering first the rising 
component, the parameters of the dipole remain essentially 
unchanged. Using now the sinking component only, the de- 
rived field parameters are very different from anything we 
obtained before. While (3 keeps its typical value (36 ±4°), 
Bp is now equal to 401 ± 26 G. This large increase in field 
strength is a direct consequence of the fact that the cool 
component of the atmosphere contributes for l/10 th only of 
the continuum flux, so that a magnetic field 10 times more 
intense than before is needed to account for the amplitude of 
Zeeman signatures. The minimum Xi reached in both cases 
(2.6 and 2.3 for the hot and cool component respectively) 
shows that no improvement in data fitting is obtained. 

As a last test, we calculate several 3D-x?-landscapes for 
different values of the vertical velocity of magnetic regions. 
The best result is obtained for a down flow of 2.3 kms" 1 (in 
agreement with the observed red-shift of the zero-crossing). 
The dipole parameters are not very affected by this velocity 
shift, and the corresponding Xr is equal to 1.4. We note that 
this velocity is quite far from the typical value of 5 kms 
we would expect in the case of a convective down-flow. We 
explore further in the next section the possibility that mag- 
netic elements may be confined in down-flows and show that 
this idea is not consistent with the basic properties of tur- 
bulent convection. 

6.2 Temperature of magnetic regions 

If the apparent red-shift of polarized profiles is linked to a 
concentration of magnetic regions in down-flows, the clas- 
sical picture of turbulent convection imposes that magnetic 
regions should be cooler than the quiet photosphere. To test 
this idea, we calculate two new sets of Stokes V LSD profiles, 
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using two sub-masks built up respectively with the blue and 
red parts of the G7 line-mask. The wavelength limit is taken 
equal to 550 nm, ensuring that both sub-masks produce line- 
profiles of similar S/N. The mean wavelengths of the blue 
and red masks are 477 nm and 587 nm respectively. In both 
line-lists, the mean Lande factor is close to 1.21. In Fig. 
we illustrate the differences observed between the averaged 
blue and red Stokes V profiles. Note that the blue profile is 
slightly expanded to compensate for a difference in equiva- 
lent width of about 7% between blue and red Stokes I LSD 
profiles. 

A difference in the amplitude of Zeeman signatures 
is visible, with peaks at 2.4 x 1CT 4 ± 2 x 10~ 5 / c and 
3.9 x 10~ 4 ± 2 x 10 _5 7 c for the red and blue signatures, re- 
spectively. If the Zeeman signatures are formed uniformly in 
the photosphere (of temperature 5550 K), no significant dif- 
ference should be observed between the blue and red profiles. 
The blue excess we can see here therefore suggests that Zee- 
man signatures are preferentially formed in hotter regions 
of the photosphere. If we exclusively link the blue excess 
to a temperature effect, the implication is that magnetic 
regions are several thousand Kelvin hotter than the quiet 
photosphere. This temperature difference is too important 
to be considered as a satisfactory explanation. We therefore 
suggest that additional (still not understood) effects must 
be taken into account to reproduce the observed blue/red 
contrast of Stokes V profiles. 

If magnetic elements producing the observed Zeeman 
signatures are hotter than the average temperature of the 
photosphere, no red-shift of the signatures should be ob- 
served, at least as long as turbulent convection is concerned. 
We therefore conclude that the apparent red-shift of Stokes 
V profiles is not related to a systematic vertical velocity of 
magnetic regions. In this context, the presence of a large- 
scale toroidal field provides a much more satisfactory inter- 
pretation of the observed shape of polarized profiles. 



6.3 Asymmetry of Stokes V profiles 

As already pointed out in Sect.|3] an asymmetry is observed 
in Stokes V profiles, with an amplitude and area of the blue 
lobe exceeding that of the red lobe. This asymmetry can also 
be readily seen in Fig. |5] where the (symmetric) derivative 
of Stokes I profile is compared to the averaged Stokes V 
profile. We note that the level of asymmetry is higher for 
the mean profile calculated with a red line-masks (Fig. |UJ. 
For the blue profile, we derive SA = 0.17 ± 0.06 and 5a = 
0.22 ± 0.07. For the red profile, SA = 0.34 ± 0.08 and 5a = 
0.35±0.09 (where SA and SA, as defined in Sect. [3] represent 
the relative asymmetry affecting respectively the area and 
amplitude of the lobes). 

In this section, we introduce an asymmetry in our syn- 
thetic polarized profiles and test whether this new parame- 
ter can affect the reconstructed field geometry. We simulate 
asymmetric local Stokes V profiles by expanding one lobe of 
the profile and reducing the other lobe by the same factor. 
The local profile is then smeared by the instrumental reso- 
lution and by the rotational broadening. The effect of this 
smearing is to decrease the amplitude of the profiles and 
generate an apparent shift of the zero-crossing wavelength 
(as described by Solanki & Stenflo 1986). The best fit of the 
mean Stokes V profile is obtained for a relative asymmetry 
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10. Same as Fig. [5] for a magnetic model including an 
dipole, a toroidal field and an asymmetry of Stokes V 



of the local profile 5a ~ SA ~ 0.11, which is then used to 
calculate magnetic models. 

Assuming a dipolar field geometry as in Sect. 14.31 the 
addition of the asymmetry factor does not lead to any sig- 
nificant impact on the dipole parameters. The Xr ls however 
significantly improved and is now equal to 1.87 (vs. 2.35 
without asymmetry). However, this relatively high \r value 
demonstrates that an asymmetry factor is not sufficient to 
"mimic" the shape of observed Stokes V profiles (Fig. 1111 . 

Assuming a field geometry composed from a global 
dipole and an axi-symmetric toroidal field, the value of B p 
keeps close to that obtained with symmetric profiles. By 
contrast, the value of B t is significantly decreased, from 
-148 ± 9 G to -112 ± 8 G. This evolution is due to the 
red-shift of the zero-crossing velocity induced by the combi- 
nation of a local asymmetry and line smearing. Here the Xr 
is also slightly improved with respect to the original model, 
with a value of 1.23. The fit of individual phases is shown in 
Fig.lTUl 



7 DISCUSSION 

7.1 Large-scale dipole 

The modeling of our data set suggests that one component 
of the large-scale magnetic field of £ Bootis A is an inclined 
dipole, with a strength of about 40 G and an inclination 
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Figure 11. Averaged observed Stokes V profile (full line) compared with mean profiles obtained from different modeling options described 
in the text (dashes). Left panels : dipolar field without and with asymmetric synthetic profiles (top and bottom, respectively). Right 
panels : dipole and azimuthal field without and with asymmetry (top and bottom, respectively). 



angle of order of 35° with respect to the rotation axis. The 
dipole parameters are only slightly varying for different basic 
assumptions in our magnetic models (Tab.|5Jl. 

Using observations collected in 1998, along with archival 
magnetic field measurements by Borra et al. (1984) and 
Hubrig et al. (1994), Plachinda & Tarasova (2000) report 
a sign reversal of the longitudinal field Bi during the stellar 
rotation, which denotes a relatively high inclination of the 
dipole with respect to the rotation axis (f3 — 87° for their 
chosen inclination i — 40° ) . It is obvious from Fig. Q that 
this behavior is not observed in our data set. A quick look 
at our series of profiles is sufficient to conclude that the sign 
of Zeeman signatures is constant over the whole rotation, 
which translates into a constant sign of Bi (Fig. 2J. Our 
own estimate of (5 ranges between 20° and 40°, depending 
on the details of the underlying magnetic model. This es- 
timate is not very sensitive to i and to the rotation period 
(Sect. |5J, which are different in our study from the values 



taken by Plachinda & Tarasova. We note that, among all 
measurements used by Plachinda & Tarasova, only six pos- 
sess a negative value in excess of la (and only one in excess 
of 2a). Five of these points come from observations gath- 
ered at the Crimean Astrophysical Observatory in 1998, and 
only one from archival measurements. One possibility is that 
negative Bi values reported by Plachinda & Tarasova come 
from intrinsic biases in their measurements. A second pos- 
sibility is that the field geometry has evolved between their 
observations and our own observing campaign. Long-term 
changes in the photospheric field would also help to inter- 
pret the apparent discrepancy between their maximum B\ 
values (about 50 to 80 G) and those reported in the present 
article (never exceeding 20 G). We come back to this last 
point in Sect. 17.31 
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7.2 Large-scale toroidal field 

A simple dipole is not able to provide a good fit to our 
time-series of Stokes V profiles. This situation is dramati- 
cally improved by the addition of an axi-symmetric ring of 
azimuthal field on the visible hemisphere of the star. A large- 
scale toroidal field is also required if we take into account 
the profile asymmetry in our magnetic model, but in this 
case the intensity of the azimuthal component is decreased 
by about 20%. 

If we do not introduce any azimuthal field component, 
the only option to obtain a similar \ 2 level is to assume that 
magnetic regions are concentrated in photospheric down- 
flows, with a vertical velocity of order of 2.3 kms -1 . We con- 
sider that this model is not adapted, for two main reasons. 
The first argument is that convective down-flows would cor- 
respond to material cooler than the quiet photosphere, while 
we deduce from the blue/red contrast of Zeeman signatures 
that the magnetic regions are, to the opposite, hotter than 
the average photospheric temperature. The second argument 
is that typical velocities of magnetic elements observed on 
the solar surface are even smaller than typical solar convec- 
tive velocities and generally do not exceed a few hundreds of 
ms _1 (Solanki 1993), which is far below the value required 
here. Since convective velocities increase with surface tem- 
perature (Gray & Toner 1985), we even expect potential 
systematic velocities of magnetic elements to be smaller on 
£ Bootis A than on the Sun, and in any case much smaller 
than 2.3 kms -1 . From these two arguments, we conclude 
that a toroidal field on £ Bootis A is the only acceptable 
interpretation of our spectropolarimetric observations. 

Azimuthal magnetic fields were previously detected on a 
small number of fast-rotating stars (e.g. Donati et al. 2003, 
Petit et al. 2004b) in striking contrast to the solar case, 
where the largest magnetic structures generally host a ver- 
tical field (except in the penumbra of sun-spots). In a few 
cases, the toroidal component reconstructed by ZDI is suffi- 
ciently structured to display a clear latitudinal dependence, 
such as rings encircling the pole at different latitudes. 

The possible detection of a toroidal field on £ Bootis A, 
obtained with a modeling technique independent from ZDI, 
suggests that large-scale azimuthal fields may be common 
on very active solar-type stars. Despite its low projected ro- 
tational velocity, the rotation period of £ Bootis A is equal 
to 6.43 d, which makes it a much faster rotator than the 
Sun. Since it has been known for a very long time that the 
efficiency of dynamo processes is correlated to the stellar ro- 
tation rate (e.g. Baliunas et al. 1995), it is not a surprise 
that we observe here magnetic characteristics already re- 
ported for intermediate rotators of rotation periods of a few 
days. 

7.3 Temporal evolution of the magnetic field 

We report a significant short-term evolution of the field ge- 
ometry during the 40 consecutive nights of data collection. 
This evolution translates into an increase of B p , /3 and B t . 
At this stage, it is not possible to decide if this changing 
geometry is due to a global evolution of the large-scale field, 
or to the continuous redistribution of active regions. Consid- 
ering the rather short timescale of the evolution, we suggest 
the second option is more likely. Such fluctuations in the 



parameters of the large-scale field are expected, since the 
"global field" we measure here is the average, over the vis- 
ible hemisphere of the star, of a potentially very complex 
and fluctuating distribution of individual magnetic regions. 
Future observations may tell us however if we observed here 
the first steps of a more global geometry change. 

One way to interpret the apparent discrepancy between 
the magnetic geometries derived in the present study and 
in the previous work of Plachinda & Tarasova (2000) is to 
invoke an intrinsic evolution of the global field. This would 
imply that the field geometry has significantly evolved in 
less than 5 yr (since 1998), with a marked decrease of the 
dipole inclination (from 87° to about 30°), together with a 
decrease in the dipole strength (from about 80 G to 40 G). 
This behavior would be very different from what is observed 
on the Sun, where the large-scale dipole is almost aligned on 
the rotation axis (at least when its strength is maximum, 
i.e. during solar minimum). Recent numerical simulations 
of magnetic flux transport processes on the solar surface 
suggest that the polarity reversal of the global dipole can 
be interpreted as a progressive 180° shift of its inclination 
(Baumann 2005). However, this short transition is occurring 
at solar maximum, when the strength of the dipole is very 
small (to a level where the dipole is probably not detectable). 
For £ Bootis A, this general picture is in contradiction with 
the apparent long-term evolution of the field geometry (sug- 
gesting that a higher inclination is associated to a stronger 
dipole). Only a monitoring of £ Bootis A over several years 
with a unique instrument, together with a field modeling 
performed with a unique procedure, will make it possible to 
determine if such non-solar field evolution is genuine. In the 
meantime, we cannot exclude that the observed differences 
with previous field measurements by Plachinda & Tarasova 
partly arise from differences in instrumental configuration 
and /or data processing. 

7.4 Nature of magnetic regions 

The large-scale field investigated in the present study results 
from an average, over the visible hemisphere of the star, of 
magnetic structures of possibly various natures. The fact 
we observe a short-term evolution of the field suggests that 
individual magnetic regions at the origin of the large-scale 
field have typical life-times of days to weeks at most. More- 
over, the blue/red contrast of Stokes V profiles indicates 
that magnetic elements we observe may be hotter than non- 
magnetic regions of the stellar surface. Both observations 
can be reconciled with what is observed on the Sun in facu- 
lae and in the network field. However, rotational variations 
in bisectors of Stokes I profiles provide a marginal evidence 
that the granulation velocity dispersion is higher for rotation 
phases at which the magnetic field is maximum. If confirmed 
by future observations (benefiting from a better spectral res- 
olution), this may imply that the general characteristics of 
magnetic regions are very different from anything observed 
on the Sun. In particular, this could suggest that the "star- 
patch" reported by Toner & Gray (1988) is associated to a 
large-scale structure of the photospheric field. 

In order to reproduce the observed asymmetry of Stokes 
V profiles, we suggest that the local polarized profiles (i.e. 
prior to instrumental and rotational smearing) must possess 
a relative asymmetry 5a w S A « 0.11. The value of 8a is 
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similar to what is observed on the Sun (e.g. Pantellini et al. 
1988). 8 A is also close to solar measurements, though slightly 
higher than typical solar values. The asymmetry of solar pro- 
files is usually modeled by the combined use of velocity and 
magnetic gradients in magnetic elements (see Solanki 1993 
for a review). The sign and level of 8a and 8 A depend on the 
nature of magnetic regions in which the Zeeman signature is 
formed (faculae and network field are considered by Solanki 
1989, sunspot penumbra by Solanki & Montavon 1993). In 
the future, a rigorous modeling of the asymmetric profiles 
of £ Bootis A may help to distinguish between the different 
possible magnetic structures at the origin of the observed 
asymmetry. 



8 CONCLUDING REMARKS AND FUTURE 
WORK 

Our direct modeling of the circularly polarized line profiles of 
£ Bootis A suggests that the large-scale field is mainly con- 
stituted of two main components, with a global dipole and 
a large-scale toroidal field. Using different magnetic models, 
we derive a dipole strength varying between 27 and 43 G, 
a dipole inclination ranging from 9° to 35° and a toroidal 
field strength spreading between 112 and 160 G. The ob- 
served range in magnetic parameters have several identified 
origins. First, the dipole parameters are affected by the addi- 
tion of the toroidal field in our models. Second, the toroidal 
field strength itself is modified by the further assumption 
that Stokes V profiles are asymmetric. Third, all field pa- 
rameters (dipole and toroidal component) are affected by a 
short-term evolution of the field geometry over the 40 nights 
of our observing window (with an increase of field strength 
and dipole inclination). 

This work confirms that azimuthal field components 
may be common on very active solar-type stars (see also 
results derived from ZDI : Donati et al. 2003, Petit et al. 
2004b). Such information cannot be extracted from model- 
ing techniques based on B t measurements, which make use 
of a small fraction only of the information available in the 
data. We demonstrated in the present study that, even for 
slow rotators with narrow line-profiles, tight constraints on 
the magnetic geometry can be obtained with a modeling of 
the full profiles. The detection of a Stokes V asymmetry is 
another example of information that can only be obtained 
using high S/N polarized profiles. 

We have engaged a long-term monitoring of £ Bootis A, 
in order to investigate the secular evolution of its global field. 
We will test in particular whether the large-scale field un- 
dergoes periodic polarity reversals (as in the solar case) or if 
fluctuations of the magnetic geometry are more erratic (as 
suggested by Baliunas et al. 1995 from a long-term monitor- 
ing of chromospheric activity). We will also follow potential 
long-term fluctuations in the blue/red contrast of Zeeman 
signatures and in the line asymmetry. 

This work opens exciting perspectives for the study of 
slowly-rotating solar-type stars with moderate to low mag- 
netic activity, for which classical Doppler imaging techniques 
are not applicable. We have demonstrated that their global 
field can be investigated with an accuracy of a few Gauss, 
and we expect even better performances with the new gen- 
eration of spectropolarimeters like ESPaDOnS/CFHT and 



NARVAL/TBL, which will cover a broader spectral domain 
and benefit from a higher spectral resolution (Donati et al. 
2005). The efficiency of multi-line techniques will be signifi- 
cantly increased with such instruments and therefore enable 
the investigation of dynamo cycles on a sample of solar-type 
stars, with the aim to estimate how various stellar parame- 
ters may impact general properties of the large-scale dynamo 
at work in active stars. 
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